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ABSTRACT 

The  problem-  of  transverse  twisting  of  a shallow  spherical 
shell  with  a small  circular  hole  is  solved,  in  generalization 
of  the  corresponding  problem  of  a flat  plate.  The  solution  is 
of  interest  as  a closed-form  solution  of  an  unsymmetrical  stress 
concentration  problem,  with  quantitative  features  depending  on 
its  boundary  layer  behavior  for  large  values  of  a relevant 
parameter.  The  problem  is  also  of  interest  as  an  example  of 
applicability  of  a previously  proposed  asymptotic  procedure  where 
interior  contributions  and  edge-zone  contributions  are  determined 
in  sequence  rather  than  simultaneously. 


On  the  Transverse  Twisting  of  Shallow  Spherical  Ring  Caps* 

by 

E.  Reissner 

Department  of  Applied  Mechanics  and  Engineering  Sciences 
University  of  California,  San  Diego 
La  Jolla,  California  92093 

Introduction.  The  original  aim  of  this  note  was  to  formulate 
a nonrotationally— symmetric  stress  concentration  problem  for  thin 
shells  which  could  be  solved  in  closed  form,  and  to  obtain  the  solution 
of  this  problem.  It  appeared  in  the  course  of  the  analysis  that  this 
stress  concentration  problem  was  also  a particularly  fitting  example 
for  the  application  of  an  asymptotic  solution  method  for  unsymmetric 
shell  problems.  Involving  the  concepts  of  interior  and  edge  zone  solu- 
tion contributions  and  of  the  concept  of  contracted  boundary  conditions 
for  the  separate  determination  of  these  contributions,  which  had  been 
proposed  sometime  earlier  [4], 

The  problem  is  as  follows.  We  consider  an  Isotropic  shallow  spherical 
shell  with  the  edges  defined  by  two  pairs  of  mutually  perpendicular  planes 
perpendicular  to  a base  plane,  with  the  corners  of  the  rectangle  in  the 
base  plane  which  is  determined  by  the  two  pairs  of  mutually  perpendicular 
planes  coinciding  with  the  corners  of  the  shell  boundary  curve.  Given 
this  configuration,  we  assume  that  the  edges  of  the  shell  are  free  of 
stress,  except  for  the  action  of  equal  and  opposite  concentrated  corner 
forces,  as  indicated  in  Figure  1.  Our  object  is  the  state  of  stress  in 
the  shell,  without  or  with  a small  concentric  circular  hole  at  the  apex. 

It  is  evident  that  a limiting  case  of  the  above  problem  is  the 
corresponding  problem  of  a flat  plate,  with  the  solution  of  the  problem 
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without  the  circular  hole  being  a special  case  of  the  problem  of 
St.  Venant  torsion  of  narrow  rectangular  cross  section  beams,  and  with 
the  solution  of  the  circular-hole  problem  being  included  in  solutions 
by  Goodier  for  a class  of  transverse  plate  flexure  problems  [2). 

In  the  present  analysis  the  plate  flexure  problem  appears  upon  • 
assuming  the  value  of  a certain  parameter  p to  be  zero.  At  the  same 
time  the  asymptotic  analysis  corresponding  to  the  procedure  described 
in  [A]  is  appropriate  for  values  of  p which  are  large  compared  to  unity. 

In  the  interim  region  of  finite  values  of  p it  is  necessary  to  obtain 
appropriate  solutions  of  the  equations  of  shell  theory,  which  In  this 
instance  may  be  taken  from  shallow-shell  theory. 

Regarding  the  physical  aspects  of  the  problem  we  find,  as  expected, 
a dominance  of  bending  stresses  over  membrane  stresses  in  the  interior 
of  the  shell  region.  On  the  other  hand,  we  also  find  that  for  suffi- 
ciently large  values  of  p we  have  membrane  stresses  in  an  edge  zone 
which  are  of  the  same  order  of  magnitude  as  the  bending  stresses  in  this 
zone,  in  such  a way  that  the  value  of  the  stress  concentration  factor 
for  this  problem  of  transverse  bending  involves  both  bending  and  mem- 
brane stresses  in  a significant  manner. 

Equations  for  Isotropic  Homogeneous  Shallow  Spher leal  Shells. 

We  consider  a shallow  spherical  shell  with  middle  surface  equation 
z * H - r*/2R,  where  R is  the  radius  of  the  shell,  H the  distance  of  the 
apex  from  the  base  plane  of  the  shell,  and  r and  6 are  polar  coordinates 
in  the  base  plane.  We  assume  that  the  shell  is  free  of  distributed 
surface  forces  and  have  then  that  tangentional  stress  resultants  N,  stress 
couples  M and  transverse  stress  resultants  0 are  expressed  as  follows 
in  terms  of  a stress  function  K and  a transverse  displacement  function 
w.  [3), 
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Use  of  appropriate  equations  of  equilibrium  and  compatibility  in 

conjunction  with  the  above,  and  in  conjunction  with  stress  strain  relations 

of  the  form  e _ - B(N  -vHa.) , etc.  leads  to  differential  equations  for 
rr  n uv 

K and  w of  the  form 


RBV*V*K  - V*w  - 0 . RDV*V*w  + V*K  - 0 
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It  is  readily  verified  that  the  solution  of  the  system  (4)  may  be 
expressed  in  terms  of  three  functions  i>  and  x In  the  form  [4], 

w - ♦ + x . K ” ♦ - RDV*x  , 


(5) 


provided  that 

V*4  - 0 . ?*♦  - 0 . V’V’x  + X-x  - 0 , (6) 

where  X*  ■ 1/R*BD. 

We  note  for  what  follows  as  expressions  for  resultants  and  couples 
in  terms  of  4,  *>  and  x 
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and  we  alao  note  the  designations  of  4 and  f as  lnextensional  bending 
and  membrane  (interior)  solution  contributions,  respectively,  and  the 
designation  of  x as  edge  zone  solution  contribution,  with  the  physical 
significance  of  the  latter  designation  depending  on  an  appropriate  relation 
between  the  length-parameter  1/X  and  an  appropriate  linear  dimension 


of  the  shell. 

The  Boundary  Value  Problem.  We  start  out  with  the  observation  that 
the  classical  solution  w - -Pxy/2(1-V)D  for  St.  Venant  twisting  of  a 
flat  rectangular  plate  as  produced  by  an  arrangement  of  concentrated 
corner  forces  P,  in  conjunction  with  an  assumption  of  no  inplane  stress, 
that  is,  in  conjunction  with  the  stipulation  K - 0.  also  satisfies  the 
differential  equations  (4)  for  shallow  spherical  shells.  Furthermore, 
this  solution  of  (4)  satisfies  the  same  corner  force  conditions  for 
a spherical  cap  with  otherwise  free  edges,  in  the  event  that  the  projec- 
tion of  these  edges  onto  the  base  plane  of  the  shell  happens  to  be 
rectangular. 
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Having  the  above  staple  solution  for  transverse  twisting  of  a 
spherical  cap,  we  ask  for  the  way  In  which  this  solution  Is  modified 
by  the  presence  of  a circular  hole  of  radius  a,  concentric  with  the 
apex  of  the  shell,  given  that  a is  small  compared  to  the  overall  dimen- 
sions of  the  cap.  Evidently,  the  boundary  conditions  for  the  free 
edge  of  this  hole  are  of  the  from 

r ’ Nrr  “ Nre  - Mrr  * Or  + r"Mre,a  " »•  <14> 


As  regards  the  boundary  conditions  along  the  outer  edges  of  the 

cap,  we  make  the  stipulation  that  for  large  r we  will  have  a homogeneous 

state  of  stress  with  cartesian  couple  and  resultant  components  - -P/2, 

jcy 

^ - Myy  - 0,  - Qy  “ - Nyy  - Nxy  - 0.  This  is  transformed,  in 

an  elementary  manner,  into  four  conditions  of  the  form 
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Closed-Fora  Solution  of  the  Boundary  Value  Problem.  The  form  of  the 
boundary  conditions  (14)  and  (15) , in  conjunction  with  the  form  of  the 
differential  equations  (4)  indicates  that  suitable  expressions  for  w and 
K will  be  product  solutions  f(r)  sin  20.  Considering  that  w and  K must 
be  as  in  (5)  and  (6),  and  deleting  at  the  outset  terms  not  compatible 
with  the  prescribed  boundary  conditions  at  infinity,  we  have  then  that 
w and  K will  be  of  the  form 
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with  four  arbitrary  constants  c , and  with  the  Kelvin  functions  ker  and 

n ? 
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kei2  subject  to  the  two  ordinary  second  order  differential  equations 

ker'2'x  + x"'kerjx  - 4x-2ker2x  - -kei2x  , (18a) 

kei'2'x  + x"'kei2x  - 4x~2kei2x  - ker2x  (18b) 

In  deriving  expressions  for  stress  resultants  and  couples  from  (16) 
and  (17),  it  will  be  convenient  to  Introduce  the  abbreviations 

ker2  - kr  , kei,  - kj  ; Ar  - x . Aa  «■  p (19) 

Therewith,  and  with  (18a, b),  we  obtain  from  equations  (1)  and  (3) 
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Introduction  of  (20)  to  (23)  into  the  boundary  conditions  (14) 
then  leads  to  the  following  set  of  four  simultaneous  equations  for  the 
determination  of  the  four  constants  of  integration  cn, 

c^tikj  - 4kt)  - c%(„k;  - 4kr)  - -6c2  , (24) 
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where  now  k^  = kj.(u),  etc. 

Upon  suitable  transformations,  this  system  of  equations  can  be 
written  in  a somewhat  simpler  form.  To  begin  with,  equations  (24)  and 
(25)  are  readily  shown  to  be  equivalent  to  the  set* 

-c,  + c,k^  - cvkr  * 0 , (24') 

2c 2 + C,ukj  - c^pk^  - 0 . (25') 

Having  (24')  and  (25'),  we  may  use  (26)  and  (27)  so  as  to  obtain 
in  place  of  these  two  equations  the  set 


-2c,  - fa  c,  + c,uk;  + c%pkj  - -1  , (26*) 

c>  - \ c.  + c,kr  + c*ki  " " ? • (27*) 


Before  evaluating  the  system  (24')  to  (27*),  it  is  useful  to  establish 
the  analytical  form  of  the  quantities  which  are  of  principal  physical 
interest.  These  quantities  are  the  edge  values  of  the  couple  Mno  and 
the  resultant  . We  obtain  a particularly  convenient  form  of  these 
expressions  by  making  use  of  equations  (1)  and  (2),  in  conjunction  with 

* Corresponding  to  the  fact  that  the  conditions  Nrr  « Nre  - 0 for  r - a 

can  be  shown  to  be  equivalent  to  conditions  K - K * - 0. 

* * 
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two  of  the  boundary  conditions  In  (14) , so  as  to  have 


M0e(a,0)  - -(l+v)DV’w(a,0)  , N00(a,0)  = V?K(a.0) 


An  introduction  of  (16)  and  (17)  into  (28)  gives,  with  the  help  of  (18a, b), 
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Having  (29a, b)  we  see,  with  the  help  of  (24')  and  (27'),  the  possibility 


of  the  further  relations 
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and  It  remains  only  to  determine  the  coefficients  c3  and  c2  from  equa- 
tions (24')  to  (27*).  We  do  this  by  first  expressing  c3  and  cH  in  terms 
of  c2,  from  (24')  and  (25'),  in  the  form 
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and  by  then  using  (26*)  and  (27')  in  order  to  obtain  the  relations 
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It  is  possible  to  simplify  the  form  of  (32a, b)  somewhat  by  making 
use  of  certain  identities  involving  Kelvin  functions  of  various  orders. 


In  this  way  we  obtain*,  upon  Introducting  (32a, b)  Into  (30a, b),  as 
expressions  for  the  significant  edge  moment  and  the  significant  edge 
resultant,  in  terms  of  zeroth  order  Kelvin  functions. 
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Stress  concentration  factors  for  bending  stresses  and  membrane 
stresses.  We  define  a bending  stress  concentration  factor  k^  as  the 
ratio  M00(a,Ti/4)/Mo  where  Mq  = M00(°°,x/4)  = -P/2.  Therewith  k^  is 
directly  given  by  the  right-hand  side  in  (33a) . 

In  order  to  obtain  the  corresponding  membrane  stress  concentration 
factor  km>  it  is  necessary  to  be  more  specific  about  the  nature  of  the 
two-dimensionally  isotropic  shell  medium.  We  shall  assume  in  what 
follows  that  the  shell  is  homogeneous  in  thickness  direction  and  have 
then  the  relation 
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We  write  further 
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and  therewith  obtain  from  (34a) 
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Stress  concentration  factors  for  small  and  for  large  values  of  p. 
Given  that  u = Aa  = a/ * Vl2(l-v*)a/ /KR,  the  limiting  case  of  a 
flat  plate  corresponds  to  the  assumption  u « 0.  We  find,  from  equations 
(33b)  and  (34b),  that  f^O)  - -1/4  and  f2(0)  ■ 0 and  therewith  from 
(33a)  and  (35c), 


4 + 4v 
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(36a, b) 


with  this  result  coinciding,  as  it  should,  with  Goodier’s  result 
for  plates,  without  consideration  of  transverse  shear  deformation  1 2 J . 

For  the  case  of  large  u,  corresponding  to  a shell  problem  with 
distinct  interior  and  edge  zone  solution  contributions  use  may  be  made 
of  appropriate  asymototic  formulas.  We  find,  by  making  use  of  certain 
known  cross-product  expansion  formulas*  that 
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Equations  (9.10.32)  to  (9.10.34)  in  [ 1 1 . 
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except  for  terms  of  relative  order  l/p?. 

Inasmuch  as  bending  and  membrane  stresses  superimpose  the  relevant 
stress  concentration  factor  for  the  most  highly  stressed  face  of  the 
shell  comes  out  to  be,  for  sufficiently  large  values  of  p, 

k - kb  + k„  = 1 + v + Vt2-  - ^ ■ <*» 


It  may  be  noted  that  the  numerical  values  of  k for  p * 0 and  for 

p = " are  not  greatly  different,  but  that  while  for  p « 0 the  stress 

concentration  is  due  entirely  to  bending,  a significant  fraction  of  it 

is,  for  1 <<  p,  due  to  membrane  rather  than  due  to  bending  action. 
Numerical  values  for  f , , f . , k.  and  k , as  a function  of  p and 

v , may  be  found  in  Table  1. 

Interior  solution  stresses  for  large  p.  The  form  of  the  expressions 
(16)  and  (17)  for  w and  k indicates  that  for  large  values  of  p the 
effect  of  the  terms  with  c,  and  c„  is  significant  in  a narrow  edge  zone 
only  and  that  outside  this  zone  the  remaining  expression  for  w is  as  if 
bending  occurred  without  stretching  and  the  remaining  expression  for  K 
is  as  if  the  state  of  stress  of  the  shell  was  a pure  membrane  state. 

We  obtain  information  on  the  state  of  stress  outside  the  narrow 
edge  zone,  and  in  particular  on  the  relative  significance  of  bending 
and  membrane  stresses,  by  determining  the  values  of  and  Mfl0  in 
accordance  with  (16) , (17)  and  the  defining  relations  (1)  and  (2) , hy 
setting  c,  - cM  = 0 in  (16)  and  (17)  and  by  then  deriving  the  relations 
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We  evaluate  (40a)  by  taking  c(  from  equation  (26) , with  c * and 
as  In  (31)  and  (32a).  Therewith  we  obtain,  except  for  terms  small 
of  higher  order 


.<i> 


= c. 


(41a) 


A corresponding  evaluation  of  (40b)  leads  to  the  relations 
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A comparison  of  (41a, b)  with  (38a, b)  shows  that  the  order  of 
magnitude  of  the  bending  stress  In  the  Interior  Is  the  same  as  the 
order  of  magnitude  of  this  stress  in  the  edge  zone,  in  such  a way  that 
the  dimensionless  edge  zone  value  1+v  decreases  to  a value  1 in  the 
Interior.  At  the  same  time  the  interior  membrane  stress  comes  out 
to  be  small  of  relative  order  1/p2  so  that,  effectively,  the  Interior 
state  of  the  shell  is  a state  of  inextensional  bending. 

Direct  Asymptotic  Solution  for  Interior  and  Edge  Zone  States. 

We  proceed  as  in  (41  to  solve  the  given  boundary  value  problem,  for 
values  of  P which  are  sufficiently  large  compared  to  unity,  through 
use  of  equations  (5)  to  (13).  Introduction  of  (7)  and  (13)  Into  the 
two  sets  of  boundary  conditions  (14)  and  (15)  then  leaves  as  conditions 
for  the  determination  of  the  two  harmonic  functions  $ and  i|i  and  of  the 
"plate  on  an  elastic  foundation"  function  x,  for  r - a. 
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\ with  equations  (42),  (43)  and  (45)  also  holding  for  r - and  with  the 
right-hand  side  of  (44)  being  replaced  by  -(P/2D)  sin  20  for  r » «>. 

We  now  note  that  when  1 <<  p we  have  the  order  of  magnitude 
relations , 


X * o(a*>r)  , x>r  = °(axirr>  . (46) 

etc.  We  use  these  for  an  asymptotic  solution  of  the  problem,  by  retaining 
in  (44)  and  (45)  the  highest  and  second  highest  order  of  magnitude  terms 
in  X,  <?2x)>r  and  VJx,  only,  that  is,  we  replace  equations  (44) 
and  (45)  by  the  abbreviated  equations 


(l-v)$ ,ir  + 7*x 


0 . 


+ (V*x)  _ 

9 * 


0 . 


(47,48) 


An  introduction  of  this  into  (42)  and  (43)  then  leaves  as  two  conditions 
for  the  determination  of  the  two  harmonic  functions  4 and 


\rr 

- ™ ^|( 

\r 

r 

-*) 

+ ij> 

♦,r0 

♦ 0 

- v*  nn 

1-v 

|/*.r 

♦ 

r 

J T T\  u 

r» 

r' 

rr 


,00 


,00 


,rr 


,0 


(49) 

(50) 


Nl  + 


N 


Note  that  upon  writing  equations  (7)  to  (13)  in  the  form  Nt 

etc  , so  as  to  distinguish  between  interior  and  edge  zone  solution  contri- 
butions, equations  (49)  and  (50)  are  equivalent  to  the  previously  derived 
contracted  boundary  conditions  for  the  determination  of  the  interior  state 

[4],  of  the  form  r?N 


rr 


R(M, 


rO 


r’N*e  - - wjr)>e  - 0. 
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Having  determined  4 and  y>,  we  subsequently  determine  the  associated 
approximation  for  the  edge  zone  function  x with  the  help  of  equations 
(48)+  , and  we  use  the  results  obtained  in  this  way  in  order  to  obtain 
from  equations  (8)  and  (12)  as  approximate  expressions  for  the  relevant 
edge  values  of  circumferential  stress  resultant  and  stress  couple 

Nee  • *,rr  + nx'x  ■ Mee  - -d-v)D*  rr  - vDv*x  , (51) 

for  r - a. 

In  order  to  carry  out  the  remaining  simple  calculations  we  write, 
consistent  with  (16)  and  (17),  in  order  to  assure  satisfaction  of  all 
conditions  at  infinity 


♦ * w1 


Pa1  sin  26  /It2.  a*  \ 

‘ -7U-v)D  l 7 iT  + c, 


* - K 


Pa*  sin  29 
2(l-v)/I5ff 


and  we  further  write 


X * 
and 


-X(r-a)/^7 


(C  cos  + C.  sin  \Zz*\ 

' /z  n ) 


v’x  - X2e-*(T-»)//7  /c  n x£la  _ 


Ic  sin  X—  - C cos  A— \ 

\ n * /i  / 


(52) 

(53) 


(54) 


(55) 


We  now  introduce  (52)  and  (53)  into  the  boundary  conditions  (49) 
and  (50)  and  obtain  as  two  equations  for  the  determination  of  c,  and  c, 

ctu2  - (1-v)  (l+6c , ) - 0 , c,ii*  - (l-v)(l-6c))  - 0 . ( 

+ We  note  that  these  equations  may  be  written,  equivalently,  as 
DV’X  * MJr  and  D<7?*>  r " r''MIo  n ‘ n-  for  r ■ a. 
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Equations  (56)  imply,  consistent  with  (32),  that 


c2nJ  - 1 - v , Cj  - 0 (57) 

Having  c2  and  c(  as  in  (57),  we  finally  obtain  and  Ck  from 
(48)  in  the  form 

C,  - -P/2DX*  , = -Cj  - P/7/DX 4  5 a * -C, 

and  therewith,  from  (51), 

»..(«•  - - rS=  • «»»(•.  i)  ' I ' 


(58) 


(59) 


The  above  expressions  for  the  edge  values  of  N00  and  M00  may  be  compared 
with  the  interior  values  of  these  same  two  quantities,  N00(a,w/4) 

■ 6P/2/DE  u*  and  M00(a,*/4)  - -P/2,  which  follow  from  (52). (53)  and  (57), 
consistent  with  the  contents  of  equation  (41) . 


Table  1 


kb 

km 

u 

f. 

*a 

v-0 

v-1/3 

v-1/2 

v-0 

v-1/3 

v-1/2  | 

0 

-0.250 

0 

1.333 

1.714 

0 

0 

IK 

.1 

-0 . 249 

0.012 

1.332 

1.599 

1.713 

0.009 

0.008 

0.007 

.3 

-0.243 

0.063 

1.321 

1.591 

1.707 

0.048 

0.041 

0.036 

.5 

-0.234 

0.122 

1.305 

1.580 

1.699 

0.092 

0.078 

0.069 

.8 

-0.219 

0.206 

1.280 

1.561 

1.684 

0.152 

0.131 

0.116 

1 

0.257 

1.263 

1.548 

1.674 

0.187 

0.162 

0.143 

2 

-0.165 

0.443 

1.198 

1.498 

1.635 

0.306 

0.271 

0.241 

3 

0.557 

1.56 

1.465 

1.609 

0.372 

0.333 

0.299 

4 

-0.114 

0.633 

1.128 

1.443 

1.591 

0.413 

0.373 

0.336 

5 

-0.098 

0.687 

1.109 

1.427 

1.578 

0.440 

0.400 

0.362 

an 

0 

1.000 

1.333 

1.500 

0.577 

0.544 

0.500 
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